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Abstract. In the radial flow of gas into a black hole (i.e. 
Bondi accretion), the infall of any entropy or vorticity 
perturbation produces acoustic waves propagating out- 
ward. The dependence of this acoustic flux on the shape 
of the perturbation is investigated in detail. This is the 
key process in the mechanism of the entropic-acoustic in- 
stability proposed by Foglizzo& Tagger (2000) to explain 
the instability of Bondi-Hoyle-Lyttleton accretion. These 
acoustic waves create new entropy and vorticity perturba- 
tions when they reach the shock, thus closing the entropic- 
acoustic cycle. With an adiabatic index 1 < 7 < 5/3, the 
linearized equations describing the perturbations of the 
Bondi flow are studied analytically and solved numeri- 
cally. The fundamental frequency of this problem is the 
cut-off frequency of acoustic refraction, below which ingo- 
ing acoustic waves are refracted out. This cut-off is signif- 
icantly smaller than the Keplerian frequency at the sonic 
radius and depends on the latitudinal number I of the per- 
turbations. When advected adiabatically inward, entropy 
and vorticity perturbations trigger acoustic waves propa- 
gating outward, with an efficiency which is highest for non 
radial perturbations 1 = 1. The outgoing acoustic flux pro- 
duced by the advection of vorticity perturbations is always 
moderate and peaks at rather low frequency. By contrast, 
the acoustic flux produced by an entropy wave is highest 
close to the refraction cut-off. It can be very large if 7 is 
close to 5/3. These results suggest that the shocked Bondi 
flow with 7 = 5/3 is strongly unstable with respect to the 
entropic-acoustic mechanism. 

Key words: Accretion, accretion disks - Hydrodynamics 
- Instabilities - Shock waves - Binaries: close - X-rays: 
stars 



1. Introduction 

The time variability of the emission from X-ray binaries 
may be related to hydrodynamic instabilities in the pro- 
cess of accretion onto a compact star. The simplest accre- 
tion flow is the stationary, spherically symmetric accretion 
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of an adiabatic gas onto a compact object, first studied by 
Bondi (1952). The stability of Bondi accretion was estab- 
lished by Garlick (1979) and Petterson et al. (1980) in the 
newtonian limit, and by Moncrief (1980) in the frame- 
work of general relativity. Kovalenko & Eremin (1998) 
suggested that non radial perturbations might reach non 
linear amplitudes in the supersonic part of the flow if the 
size of the accretor is small enough. This effect was invis- 
ible in the 3-D numerical simulations of the Bondi flow 
by Ruffert (1994a), which confirmed its stability. By con- 
trast, a strong instability appeared in the numerical sim- 
ulations of Bondi-Hoyle-Lyttleton accretion when a bow 
shock is present (e.g. Ruffert 1994b for 7 = 5/3). The 
Bondi-Hoyle-Lyttleton accretion is named after the pio- 
neering works of Hoyle & Lyttleton (1939) and Bondi & 
Hoyle (1944) concerning an accretor moving at supersonic 
speed with respect to the gas. Foglizzo & Tagger (2000) 
(hereafter FT2000) suggested that this instability might 
be explained by what they called the entropic-acoustic in- 
stability. This generic mechanism of instability in shocked 
converging flows is based on the cycle of entropic and 
acoustic waves in the subsonic region of the flow between 
the shock and the sonic surface surrounding the accretor. 
The advection of entropy perturbations towards the ac- 
cretor produces outgoing acoustic waves which propagate 
towards the shock. Perturbed by these acoustic waves, 
the shock produces new entropy perturbations, thus clos- 
ing the entropic-acoustic cycle. According to FT2000, this 
cycle is unstable (i.e. the amplitude of the new entropy 
perturbation exceeds the initial one) if the sound speed 
at the sonic point is much larger than the sound speed 
at the shock. This mechanism might play a destabilizing 
role in various astrophysical environments such as wind 
accretion or disc accretion. In this series of two papers, 
we describe the entropic-acoustic instability in the Bondi 
flow, including a spherical shock as an outer boundary 
(as in Foglizzo & Ruffert 1997). The study of this rather 
academic configuration serves two purposes: 

(i) it is the first detailed description of the entropic- 
acoustic instability in a specific flow, beyond the general 
concepts and approximations of FT2000. The Bondi flow 
is simple enough to allow for an analytical treatment of 
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3D perturbations, including non radial perturbations of 
entropy and perturbations of vorticity. 

(ii) understanding the instability of the Bondi flow 
with a spherical shock might be the first step towards 
a comprehension of the instability of non radial shocked 
accretion flows such as supersonic Bondi-Hoyle-Lyttleton 
accretion. 

In this first paper, we investigate the key part of the 
entropic-acoustic instability: the excitation of outgoing 
acoustic waves by the advection of entropy and vorticity 
perturbations. As explained in FT2000, the efficiency of 
this process is directly linked to the efficiency of the refrac- 
tion of acoustic waves, which is also studied. The efficien- 
cies of these processes depend a priori on the frequency u> 
of the perturbation, its wavenumbers in spherical coordi- 
nates (l,rn), and the adiabatic index 1 < 7 < 5/3 of the 
gas. The accretor is a compact object such as a black hole, 
absorbing all the matter falling inside the sonic radius. 
The shock, which plays no role in this first paper, is in- 
troduced in the second paper of this series (Foglizzo 2001) 
in order to close the entropic-acoustic cycle and enable a 
description of the eigenmodes of the global instability. 
The present paper is organized as follows. Linearized equa- 
tions are established in Sect. 2. The refraction of acous- 
tic waves is studied in Sect. 3. The excitation of acoustic 
waves by the advection of entropy and vorticity perturba- 
tions is examined in Sect. 4. 



2. Linear perturbations of pressure, entropy and 
vorticity in the Bondi flow 

2.1. Properties of the unperturbed flow 

We consider the radial accretion flow of a gas with con- 
stant sound velocity and density at infinity, acceler- 
ated towards a point-like accretor of mass M, which is to- 
tally absorbing. The equations describing the unperturbed 
Bondi flow are recalled in Appendix A. In what follows, 
densities and velocities are normalized by Cqo and , and 
distances are normalized to the Bondi radius GM/c 2 ^. Af- 
ter this normalization, the unperturbed flow depends on 
a single parameter 7, the adiabatic index of the gas. The 
sound velocity c s at the sonic radius r s in the Bondi flow 
diverges for 7 = 5/3: 
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The "natural" sonic frequency c s /r s is directly related to 
the Keplerian frequency luk at the sonic radius: 
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2.2. Differential equation for small perturbations 

The Euler equation, the equations of conservation of en- 
tropy and mass in the Bondi flow are linearized for small 
perturbations in Appendix B. After a Fourier transform 
in time, the perturbations are projected onto a basis of 
spherical harmonics Y™, where l,m are the numbers as- 
sociated to the latitudinal (9) and longitudinal (ip) angles. 
Kovalenko & Eremin (1998) studied the same set of equa- 
tions in order to analyze the asymptotic behaviour of their 
solutions near the origin r — 0. By contrast, we are inter- 
ested in the effect of linear perturbations on the subsonic 
part of the flow (r > r s ). A single differential equation of 
second order is obtained in Appendix B.2, with the fol- 
lowing structure: 



ai- 



d Sp 



a - 



Sp 



= b SS R + bi5K R . 



(4) 



&_5p 

dr 2 p ' 1 dr p p 

The coefficients a,Q,ai,bo,bi are functions of r which de- 
pend on the flow velocity and sound speed radial profiles 
in the unperturbed Bondi flow. Their rather lengthy ex- 
pressions are given in Appendix B, Eqs. (B.21 to B.24). 



The homogeneous equation corresponding to the left hand 
side of Eq. (|4|) describes the propagation of acoustic waves 
in the Bondi flow and their refraction by the sound speed 
gradients. The important feature of Eq. ([!]) is the presence 
of source terms on the right hand side: 

(i) SS is the perturbation of entropy. For the sake of 
simplicity, the ratio of the molecular weight \i to the gas 
constant 1Z is set to n/TZ = I throughout this paper, with 
no loss of generality. 

(ii) SK involves the curl of the vorticity vector Sw = 
V x Sv projected along the flow velocity: 



SK = r 2 v ■ (V x Sw) +/(/ + !)( 



,5S_ 

7 



(5) 



The radial dependence of the three components of the vor- 
ticity vect or Sw i s explicitely integrated in Appendix B.l 
(Eqs. B.5 to B.7). Nevertheless, the vorticity contributes 
to the excitation of acoustic waves only through the ra- 
dial part of its curl, as it appears in Eq. (|^). In an adia- 
batic flow, bot h SS and SK are conserved when advected 
(see Appendix B.f). A more compact mathematical for 



mulation is obtained in Appendix B.2. 2 by writing the 
differential equation satisfied by the perturbation / of the 
Bernoulli constant. 
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where the frequency uii is defined as: 
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Fig. 1. Refraction coefficient \1Z S \ 2 for 7 = 1.33, I = 
0,1,2,3. The dotted line corresponds to the leading be- 
haviour of \H S \ 2 for lo > ujf ut (Eq. |l|) 



The structure of this differential equation is the same as 
Eq. (0), with much simpler coefficients. The right hand 
side of Eq. (||) contains two source terms proportional to 
8Sr and SKr. The pressure perturbation is related to / 



by the following equation derived in Appendix B.l 



1 dSp 
~p~dt 



D/ 
Dt ' 



(8) 



3. Acoustic waves in the Bondi flow: the 
refraction cut-off 

3.1. Acoustic refraction and regularity at the sonic point 

The homogeneous equation associated with Eq. (](]), de- 
scribing the propagation of acoustic waves is: 



,d_ f l-M 2 d 
dr V v dr 



\-M 7 



JR. 1-M 2 " 



0(9) 



This equation is of second order, which reflects the fact 
that pressure perturbations can be decomposed into ingo- 
ing and outgoing acoustic waves. Such an identification is 
possible only in the region of the flow where the length- 
scale of the flow gradients is longer than the wavelength 
of the perturbation, where a WKB treatment of the wave 
propagation can be used (see Appendix C). This is always 
the case far from the accretor in the Bondi flow, since the 
density and sound speed are uniform at infinity. We may 
define 8p + and 8p~ as the solutions of the homogeneous 
equation propagating respectively inward and outward at 
infinity. Following FT2000, the index + (resp. — ) refers to 
waves propagating in the same (resp. opposite) direction 
as the flow. The acoustic flux carried by these acoustic 
waves is defined by Eq. (6) of FT2000: 
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Let us normalize Sp^ so that they carry the same acous- 
tic flux at infinity. Any linear combination Sp of the two 



independent solutions Sp^ , Sp , with complex coefficients 
a, (3, is also a solution of Eq. (||): 



Sp = aSp + (3Sp^ 



(11) 



There is, however, only one ratio a/f3 such that Sp is regu- 
lar at the sonic radius. Indeed, a local analysis of the solu- 
tions of Eq. @ for 7 < 5/3, at the sonic point M. — 1 us- 
ing Frobenius series, reveals that there is only one solution 
crossing the sonic point regularly, the other exhibiting a 
logarithmic singularity (Bender & Orszag 1978, Chap. 3.3, 
p. 68). This regularity condition reflects the fact that an 
incoming acoustic wave PSp^f is only partially transmit- 
ted to the accretor: an outgoing acoustic wave of ampli- 
tude aSpy is refracted out. Using the same notations as 
FT2000, the fraction \TZ S \ 2 of refracted acoustic flux is 
deduced from Eqs. dlOl) and (R~T|) : 



l^l 2 = 



(12) 
(13) 



The refraction coefficient is therefore directly related to 
the regularity of the solution at the sonic point. 

3.2. Numerical calculation of the homogeneous solution 

The differential equation (^|) is integrated using a Runge- 
Kutta implicit method, from the sonic point to the WKB 
region. The refraction coefficient is then computed by 
identifying the ingoing and outgoing waves in the WKB 
region. If 7 < 5/3, a Frobenius expansion is used in the 
vicinity of r s to start the integration away from the singu- 
larity. The calculation for 7 = 5/3 and I = is similar to 
the case 7 < 5/3, except that the sonic point is at r = 0. If 
7 = 5/3 and I > 1, the homogeneous solution presents an 
essential singularity at the sonic point, where one branch 
converges to zero, and the other diverges. Nevertheless, 
an asymptotic expansion of the converging solution (see 
Eqs. C.32j and |C.33 ) allows a numerical integration from 
the vicinity of the origin towards infinity. As expected 
by Kovalenko & Eremin (1998), no overreflection occurs 
(\TZ S \ < 1). The typical frequency dependence of \7Z S \ 2 is 
shown in Fig. |l| for 7 = 1.33. 

3.3. Cut-off frequency for the refraction of acoustic waves 
for 7 < 5/3 

If 7 < 5/3, the accretor is surrounded by a supersonic re- 
gion from which acoustic waves cannot escape. Only those 
waves with a high enough frequency may penetrate this 
region, because lower frequency waves are refracted out 
before the sonic radius. This leads us to expect the refrac- 
tion \R S \ 2 to decrease to zero at high frequency, and be 
close to unity at low frequency (see Fig. [l]). The refraction 
cut-off u>f ut of a wave of latitudinal order I is defined as 
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ca 1 (r) = 2 1 ' 2 (c 2 -v 2 ) 1/2 /r 




0.1 r sm 1 r ( ra ) 10 100 

Fig. 2. Typical radial dependence of u>\ (here 7 = 1.33), 
showing the turning point rt of a non radial acoustic wave 
I = 1 below the cut-off frequency u> < wj ut 



the intermediate frequency corresponding to a refraction 
of half of the incoming acoustic flux: 



\n. 



(14) 



Besscl functions are used in Appendix D to approximate 
the homogeneous solution at high frequency, in order to 
obtain the leading order of the asymptotic behaviour of 
\H S \. This latter decreases exponentially above ujf ut (see 
FigJ): 



0J o 



c s \M B 



(5 — 37)7r 



(15) 
(16) 



3.3.1. Non radial acoustic waves I > 1 



The turning point of non radial perturbations appears ex- 
plicitly on Eq. (^|) . It corresponds to the radius r t solution 
of u = wi(rt), where u>i is defined by Eq. (^) (see Fig. |J). 
For 7 < 5/3, this turning point is double when the fre- 
quency equals the maximum of the function uji(r t ). The 
refraction coefficient is then equal to 1/2, according to 
Bender & Orszag (1978, Chap. 10.6 p. 524): 



J i>i 



Max {l .> rs} {uji} , 

i(i + iy 



(17) 
(18) 
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Fig. 3. The cut-off frequencies Wg ut and u;J ut correspond 
to \K S \ 2 = 0.5 for I = and I = 1 respectively. They 
cannot be distinguished from their fittings (Eqs. [jl] and 
pl| ) on this plot. These frequencies are compared to the 
sonic frequency c s /r s , the Keplerian frequency at the sonic 
radius u>k, and Wqo. Frequencies are expressed in units of 

(or 7, or r s ) 



4/GM, as functions of c s /c, 



of 



J i>\ 



on the adiabatic index 7. 



u)™\ should scale as 



l x / 2 {l + l) 1 / 2 /(5 — 37) 5 / 4 according to the asymptotic be- 
haviour of (1 — M 2 ) 1 ^ 2 c/r oc CsiMs/rs) 1 / 2 when 7 — > 5/3, 
deduced from Eqs. (A. 12 to A. 14). A refined comparison 



with the results of numerical calculations of W; c >* leads to 



ufZ\~i*(i + iy- 



37 



(19) 



over the range 1.4 < c s /coo < 100 considered numerically. 
The radius r% such that wj(ri) = tof ut is the minimum 
value of the turning radius r t (see Fig. ||). ri(^) is inde- 
pendent of I, m and lies just outside the sonic surface, and 
is well fitted by r% ~ V2r s for 7 close to 5/3. 

3.3.2. Radial acoustic waves I = 



, .cut 
UJ Q 



is computed numerically in Fig. |^, and corresponds 
to a lower frequency than wj ut . The value of w§ ut can be 
estimated analytically for 7 close to 5/3 from the criterion 
of validity of the WKB approximation of acoustic waves, 
derived in Appendix C: 



to > (1 -M 2 )c 



dlogv 2 



Or 



(20) 



This result was checked by measuring numerically the re- Since the acoustic flux is conserved in the region of valid- 
fraction of acoustic waves. Fig. || shows the dependence ity of the WKB approximation, any refraction of acoustic 
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waves must occur outside this region. We deduce that w§ ut 
scales as the maximum of the function on the right hand 
side of Eq. (pp[). When 7 is close to 5/3, this maximum 
scales as (1 — M 2 )c/r cx — M s c s oc (5 — 37) -1 . This is 
confirmed by numerical calculations of w§ ut which indicate 
that in the range of sound velocities 1.4 < c s /coo < 100 
studied numerically, the ratio WQ Ut /woo is constant within 
the numerical errors: 

<j™*~ °' 85 (21) 
W ° ~5 = 37 (21) 

The scaling of w§ ut should be compared to the two other 
particular frequencies of the flow near the sonic point, 
in Eq. ([l9|) and the local Keplerian frequency luk- 



(5-37)1 



(22) 



and thus wg ut <C w£ ut <C u>k f° r 7 close to 5/3. The values 



of Wn 



computed numerically are compared to c s /r s 



and to the Keplerian frequency in Fig. |^. 



3.4- Refraction of acoustic waves for 7 = 5/3 

In a Bondi flow with 7 = 5/3, acoustic waves are perfectly 
refracted at any frequency: 

|£ s |=lfor 7 =! (23) 

This can be deduced from the limit of flows with 7 close 
to 5/3, because the cut-off frequency cuf ut tends to infin- 
ity when 7 — > 5/3. This is not surprising, since an ingoing 
acoustic wave of any frequency always ultimately meets 
"a wall" , i.e. a subsonic region where the scale of the flow 
gradients are much shorter than the wavelength. 
The WKB analysis at high frequency (Appendix C) indi- 
cates that ingoing and outgoing acoustic waves are cou- 
pled together at ro oc 1/uj for I = 0, and at their turning 
point r"t oc 1/ojs if I > 1. The continuity between the flow 
7 = 5/3 and flows with 7 close to 5/3 is checked by re- 
marking that ro ~ r s for ui cx (5 — 37) _1 and rt ~ r s for 
oj cx (5 — 37)~ 5 / 4 . Comparing this to Luf ut in Eqs. ( pl"l) 
and (|l9|), the cut-off frequency is simply interpreted as 
the maximum frequency such that the coupling of acous- 
tic perturbations occurs in the subsonic flow. 

4. Acoustic efficiency of entropy and vorticity 
perturbations 

4-1. Analytical formulation 

According to Eq. (^), it is natural to define two dimen- 
sionless complex acoustic efficiencies Qs, Qk associated 
respectively to the two source terms SSr and 5Kr. These 
definitions are made so that the acoustic flux F~ of an 
outgoing sound wave triggered by the perturbations of en- 
tropy and vorticity, measured far from the accretor is: 

2 



F~ 



Mod 



QsSS R +Q K - 



5K f 



(24) 
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Fig. 4. Acoustic efficiencies \Qs\ 2 and |Qif| 2 depending 
on the frequency for 7 = 1.33, I — 0, 1, 2, 3 
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Fig. 5. The frequencies which maximize \Qs\ are dis- 
played in units of the refraction cut-off w ; cut , for I 
0, 1, 2, 3 as functions of c s /c D 



The integral expressions for Qs and Qk are derived in 
Appendix D: 



27W Jr B 
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v{l-M 2 ) 2 d 
M 2 dr 



foe 



J 1-Ai 2 v 



Qk 



% r f Q iu rii^: 



2^ 



dr. 



dr, (25) 
(26) 



It is remarkable that both acoustic efficiencies are indepen- 
dent of the longitudinal number m. Qs and Qk depend 
only on uj and /. 

4-2. Numerical calculations 

The step of the numerical integration must be able to 
resolve both acoustic and entropy waves. Far from the 
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Fig. 6. Acoustic efficiency | Qs | 2 of entropy perturbations Fig. 7. Acoustic efficiency |Qk| 2 of vorticity perturba- 

for 7 = 5/3, I = 0, 1,2,3. The full line corresponds to the tions for 7 = 5/3, / = 1,2,3. The full line corresponds 

numerical calculation, the dotted line corresponds to the to the numerical calculation, the dotted line corresponds 

asymptotic estimate at high frequency (Eqs. and k9) to the asymptotic estimate at high frequency (Eq. p0|) 



accretor, the wavelength of the acoustic perturbation is 
asymptotically constant (~ 271000/0;), whereas the wave- 
length of the entropy perturbation decreases to zero as 
2ttv/uj oc 1/r 2 . It is therefore useful to transform the inte- 
grals in Eqs. (^5|) and (^6|) to accelerate their convergence. 
This is done in Appendix D by successive integrations by 

5/3, 



If 7 



parts, thus obtaining Eqs. ( p.6| ) and (D.19). 
the integrals are transformed near the origin in order to 
correctly treat the singularity. 

4-3. Asymptotic expansion of \Qs\ and \Qk\ at high 
frequency for 7 = 5/3 

Special properties of the Bondi flow 7 = 5/3 are related 
to the position of the sonic radius, at the central singular- 
ity. In particular, both cut-off frequencies wg ut and w£ ut 
are infinite, which means that at any frequency, incoming 
sound waves are totally refracted outward: 



1. 



(27) 



\Qs\ is estimated in Appendix E.l at high frequency for 
radial perturbations: 



|fis|i=o~f(| 



(28) 



where the Gamma function satisfies T(2/3) ~ 1.354. 
The case of non radial perturbations is derived in Ap- 



pendix E.2, where we define L 2 = 1(1 + 1): 



\Qs 



\Q 



K\1>1 



2 / loL 2 
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3 { uL- 





(2L 












(2L 
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(29) 
(30) 



The asymptotic expansions in Eqs. (|28|), (|29j) and (|30|), 
which involve no free scaling parameter, are successfully 
compared to numerical calculations of | Qs | and | Qk | in 
Figs. ^ and [?]. The radial case of Eq. (|2^) is recovered by 
taking the continuous limit / — > in Eq. (|29|). Note that 
the asymptotic behaviour of the Bessel function K2 for 
large arguments (L ^> 1), 



3?r 
4X 



(31) 



makes the acoustic efficiencies decrease rapidly for high 
order I. Remembering that the frequency of the entropy 
perturbation is directly related to its radial wavelength 
(A ~ 2itv/uj), our calculation of \Qs\{uj,1) indicates that 
the sound produced by the advection of an object of given 
mass in the Bondi flow depends a lot on its shape. Virtu- 
ally no sound comes out if I is large and u> is small (e.g. 
a stretched spaghetti). By contrast, the shape producing 
the highest acoustic flux corresponds to a small I and a 
large uj (e.g. a wide thin saucer). 

At high frequency, the coupling of entropy and vorticity 
perturbations to acoustic waves comes essentially from the 
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Fig. 8. Maximum acoustic efficiency \Qs\ 2 as a function 
of Cg/coo (i.e. of 7 < 5/3) for I = 0, 1, 2, 3. The dotted lines 



are deduced from the curves \Qs\ 7 
where lu is replaced by w ; max ~ 0.8w ; cut 
sponding to I = are marked with circles 



5/3 (Fig. |), 
Curves corre- 



region r e g oc 1/oj 2 / 3 (i.e. the main contribution to the in- 
tegrals in Eqs. and ^(| computed in Appendix E). The 
radial scaling at high frequency is therefore ro <C ft <C r e g . 
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Fig. 9. Acoustic efficiencies |Qs| 2 an d \Qk | 2 depending on 
the frequency for 7 = 1.665 (i.e. Cs/c^ = 20), I = 0, 1,2. 
Frequencies are expressed in units of c^/GM. These ef- 
ficiencies are very close to those obtained for 7 = 5/3 
(dotted lines), truncated after the cut-off frequency u; ; cut . 



for / = 1 than for I = (Eqs. [TJ and [21]). This is illus- 
trated in Fig. ^, where the acoustic efficiencies computed 
for Cg/coo = 20 (7 = 1.665) are very close to their values 
for 7 = 5/3, up to the refraction cut-off. In this figure, the 
efficiency \Qs\ is higher for I = than for I = 1 in the 
limited range of frequencies 30 < w < 120, but the overall 
maximum is reached near wj ut ~ 2300 for 1 = 1. 



4-4- Case of accretion flows with 7 close to 5/3 

For 7 close to 5/3 and at frequencies lower than the cut off 
frequency, we can use the approximation of \Qs\ obtained 
for 7 = 5/3. According to Fig. |], the maximum of \Qs\ is 
reached at the frequency w ; max , with 0.5 < w" lax / 'ujf ut < 1. 
The asymptotic behaviour of the maximum value of \Qs\ 
is therefore estimated by extrapolating Eqs. (^|) and ( p9[ ) 
to the cut-off frequency: 

MaxIQ^o-i^^g), (32) 
oc r^K (33) 
Max|e*|i>i ~ ^K ut )*^| (jf) , (34) 

oc Ltf s 0^rr*. (35) 

It is interesting to note in Fig. || that the most efficient 
entropy perturbations at a given frequency when 7 < 5/3 
are non radial (I = 1), whereas radial ones are more effi- 
cient for 7 = 5/3 (Fig. ^). This is directly related to the 
fact that the critical refraction frequency is much higher 



5. Conclusions 

5.1. Summary 

The acoustic response of the Bondi flow to pressure, en- 
tropy and vorticity perturbations has been studied thor- 
oughly. Let us summarize the main results: 

(i) the critical refraction frequency for acoustic waves 
(jf ut is quite different from the "natural" frequency c s /r s 
one might think of at first glance. Moreover uif ut is also 
very different for radial and non radial perturbations. The 
two fundamental frequencies WQ Ut , uf ut of this flow near 
the sonic point were accurately determined (Eqs. [l9] and 
f2l|) , and are ordered as follows: wg ut < wj ut < ujk, with 

oT oc Max{(l-M 2 )-}oc (^-) C ~, (36) 
I r J \Ologr J s r s 

1 

W J ut oc Max|(l-X 2 )^-}oc -■ (37) 

I ri \d\ogr J s r s 

(ii) two quantities conserved in the linear approxima- 
tion of the perturbed Bondi flow are source terms for 
acoustic waves: the entropy perturbation 8S, and a quan- 
tity 5K related to both vorticity and entropy perturba- 
tions defined by Eq. (||). 
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(iii) the acoustic efficiencies Qs and Qk of entropy and 
vorticity perturbations are independent of the longitudinal 
number m. 

(iv) \Qs\ and \ Qk\ are highest for I = 1 perturbations 
if 7 < 5/3. 

(v) \Qk\ is highest at relatively low frequency, and is 
always rather moderate (|Qif| < 1). 

(vi) \Qs\ is highest at high frequency near the refrac- 
tion cut-off, and can be much larger than one if 7 is close 
to 5/3, i.e. if c s /coo ^> 1. This confirms the physical ar- 
gument of FT2000 which stressed the importance of the 
enthalpy increase between the outer part of the flow and 
the sonic point. 

(vii) asymptotic scalings of | QsU \ Qk\ were obtained 
at high frequency for 7 = 5/3 (Eqs. |28| to |3fj|), in excellent 
agreement with numerical calculations. 

The analytical expressions obtained for a;™* (7), and 
for \Qs\ and \Qk\ at high frequency for 7 = 5/3, are 
accurate enough (see Figs. ||, || and 0) to be used as 
tests of the accuracy of 3D hydrodynamical codes used 
for Bondi accretion, and particularly their limitations at 
high frequency. 

These results suggest that if a shock were present as 
outer boundary, the entropic-acoustic cycle would be un- 
stable with respect to high frequency waves if the sound 
speed at the sonic point is high enough, i.e. if 7 is close 
to 5/3. Indeed, an entropy perturbation 5 Si advected 
towards the accretor triggers an outgoing acoustic flux 
F~ = M c 2 \Qs\ 2 \SS\ 2 . These acoustic waves propagate 
against the flow until they reach the shock, where they 
produce new entropy perturbations 5S2 (see FT2000), 
with SS2 > SSi since \Qs\ ^S> 1. The eigenmodes corre- 
sponding to this unstable cycle will be investigated in the 
second paper of this series (Foglizzo 2001). The fact that 
the most unstable simulations of Bondi-Hoyle-Lyttleton 
accretion also correspond to 7 = 5/3 (Ruffert & Arnett 
1994, Ruffert 1994b) already encourages us to look for 
an extrapolation of these results to the case of non-radial 
shocked accretion flows. 

5.2. Discussion 

5.2.1. Upper bound for |Qs| 2 imposed by the 
Schwarzschild radius 

In the flow 7 = 5/3, the fact that the outgoing acoustic 
flux produced by entropy perturbations diverges at high 
frequency (Eqs. |28| and p9| ) is rather surprising, and di- 
rectly related to the hypothesis of a point-like newtonian 
accretor (r s = 0). In this section we show that | Qs| 2 is nat- 
urally bounded if the finite size of the accretor is taken into 
account. Simple arguments are used in order to roughly 
estimate this upper bound in the case of an accreting black 
hole, without computing the exact relativistic corrections 
near its horizon. 



Let us denote by rs c h = 2GAf/cf ight the Schwarzschild ra- 
dius of the black hole, where cn g ht is the speed of light. 
With a sound speed at infinity Coo <§; cn g ht, Eq. (Q) guar- 
antees that relativistic effects are negligible except for 7 
close to 5/3. Indeed Petrich et al. (1989) checked numeri- 
cally that relativistic corrections are small for 7 = 1.1 and 
7 = 4/3. If 7 = 5/3, the sonic point is bound to lie at a 
few Schwarzschild radii: 



V 



> 1. 



^Sch 



(38) 



Our description of non radial acoustic waves in the flow 
7 = 5/3 is not affected by relativistic corrections as long as 
their turning point r t is far enough from the horizon. Since 
r t ~ (L/w) 4 / 5 at high frequency (Eq. |C.34 ), a natural 
relativistic cut-off is introduced for r t ~ r/ro,^: 



0.25 



1(1 + 1) 



Clight 



GM 



(39) 



The efficiency |Qs| 2 of the generation of non radial acous- 
tic waves close to this cut-off frequency is deduced from 
Eq. <M: 



IS* 



or 

5" 



2 



3 



Clieht 



(40) 



The numerical values of the cut-off frequency and maxi- 
mum acoustic efficiency of the entropy perturbations I = 1 
are thus typically: 



-^o I max 5 



100 

qi 
10 4 



2 ( M 



lOOkm/s/ \10AL 



Hz, 



100 km/s 



(41) 



(42) 



Therefore relativistic effects set an upper bound on the ap- 
parent divergence of \Qs\ at high frequency for 7 = 5/3. 
Eq. ( [42| ) indicates that entropy perturbations may excite 
low degree acoustic waves very efficiently, up to non lin- 
ear amplitudes if the component I = 1 of the entropy 
perturbations exceeds a few percent. This efficiency is 
much smaller for inhomogeneities with a small angular 
scale, since \Qs\ decreases exponentially when / increases 
(Eq.©. 

5.2.2. Non adiabatic limitations 

Non adiabatic processes are discussed by Chang & Os- 
triker (1985) in the context of Bondi accretion. The effi- 
ciency \Qs\ would be strongly affected by a strong ther- 
mal conduction which smoothes out entropy perturba- 
tions. This effect would be more pronounced for short 
wavelength perturbations, which are also the most effi- 
cient for 7 close to 5/3. Without going into the details 
of such processes, our understanding of the acoustic effi- 
ciency allows the following remark: if the inner region of 
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the flow were dominated by non adiabatic processes, the 
optimal frequency for \Qs\ would be at most reduced to 
the maximum frequency of acoustic waves refracted in- 
side the adiabatic region of the flow. Even when there is 
no shock, the mechanism of excitation of acoustic waves 
from the advection of entropy perturbations opens inter- 
esting perspectives to the problem of time delay in Cygnus 
X-l, approached numerically by Manmoto et al. (1996). 
Interpreting their results in terms of acoustic efficiency 
of entropy perturbations suggest that this physical pro- 
cess survives the inclusion of rotation, radiative cooling 
and viscous heating. The study of effects deserve a careful 
analysis which is beyond the scope of the present paper. 



r,v,c can be expressed as explicit functions of the variable x 

(A.10) 
(A.ll) 
(A.12) 
(A.13) 
(A.14) 



defined 


by: 


x = 


1 

1" ' 


dr 


1 - M 2 


dx 


2 ' 




(x - l) 2 (9x 4 + 2x 2 + 1) 


r = 


6a; 3 (3a; 2 - 2x + 1) 


2 

V = 


3 3x 2 - 2x + 1 


(x - l) 2 9a: 4 + 2a; 2 + 1' 


2 

C = 


3a; 4 3x 2 - 2x + 1 


(x - l) 2 9a: 4 + 2a; 2 + 1' 
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Appendix B: Linearized equations for 
perturbations 

B.l. Second order differential system 

The Euler equation is written as follows: 



Appendix A: Description of the unperturbed 
Bondi flow 

A.l. 7 < 5/3 

The Bernoulli equation and the conservation of mass are writ- 
ten after normalizing velocities to the sound speed c x at in- 
finity, and distances to the Bondi radius GM/c 2 ^ 



dv „ /V c 2 GM 

— + WXV + V [ — + 

at V 2 7 — 1 r 



2 V7 S 

= c V — . 

7 



(B.l) 



Projecting the Euler equation onto the flow velocity, we obtain 
the equation of evolution of the Bernoulli constant: 



d _\ v 2 c 2 



_ GM \ _ ldp 
r J p dt' 



(B.2) 



-r oc -1 



7 + 1 
1 \ 2(7-1) 



(I) 



(A.l) 
(A.2) 



5-37 

By taking the radial derivative of these equation, we obtain 



M 2 
1 



2 - 



<91ogc _ 7 — I M 2 

d log r 2 1 

d log v __ 1 

dlogr ~ 1 - M 2 



At the sonic point, r s ,c s are define d by Eqs. (|lj) and 
Taylor expansion of Eqs. (A.2) and (A.4) leads to: 



diogv 


2 


dlogr [rs) 


7+1 


dlogc 


7-1 


d log r s 


7+1 






dlogr 





7 - 1 + 



5 — 37 \ 2 



5 — 37 \ 2 



(A.3) 
(A.4) 

D- A 

(A.5) 
(A.6) 
(A.7) 



A.2. 7 = 5/3 



In accretion flows where 7 = 5/3, Eqs. (A.l) and (A.2) become: 

..2 Q 



u i^/2 

— C 

2 2 K 



r 



2 3 
r vc 



(A.8) 
(A.9) 



By combining the curl of Euler equation and the mass conser- 
vation, we obtain the equation for the evolution of the vorticity 
w: 



h (v ■ V)— = — • V )v + -Vc x V-. 

dt p P \P I P 7 



(B.3) 



We make a Fourier transform in time of the linearized equations 
describing the radial flow perturbed in spherical coordinated 
(r,6,<p). The linearized equation of entropy conservation can 
be directly integrated: 



SS = 8S R e 



Jr " . 



(B.4) 



The vorticity equation (B.5) can also be integrated when lin- 
earized: 



w r = 



we = — 

rv 

1 

w v = — 
rv 



(R\ 2 . , tof*: 
( — I {w r ) R e Jr " 

RvR(wg)n 



2 - c 2 o ss R 



sin 9 dip 7 

r, I \ 1 / 2 2 N 9 SSr 

Rvr{w v ) r + (c -c r ) — 

OO 7 



Jr « , 
; w f r ir. 

JR " . 



(B.5) 
(B.6) 

(B.7) 



Thus the product r 2 w r is conserved during advection, as re- 
marked by Kovalenko & Eremin (1998). In order to write the 
linearized Euler equation in the simplest form, let us define the 
two functions /, g as follows: 



/ = v Sv r + 
Sv, 



+ 



7-1 
2 5c 
7 — 1 c 



: Sc. 



(B.8) 
(B.9) 
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/ is the perturbation of the Bernoulli constant, which is d i 
rectly related to the pressure variations according to Eq. ( |B.2j ) 
g is related to the perturbation of the mass accretion rate and Eq 
entropy through: 



Using this equation in Eqs. (B.8- B.9 j ) en ables us to write (/, g) 
as functions of Sv r and 5p in Eqs. flB.ls| -B.19), and thus obtain 
with: 



9 = 



SM 
M 



SS. 



(B.10) 



a log ( c 2 



V 



2iu>v 



dr 



By linearization of Eq. (B.l), we express the velocity compo- 
nents (8vg,dVp) as follows: 



vA 

2 



L —— - --- 



1 d ( c 4 dM 2 



Svg 
V 

SVtn 



^ + 1 9 f 

iuj iuirv 88 



lui ilotv sin f 



+ • 



ilotv ' 
' d_ 



d SSr iu: 
-e 



J R v 



7 

2 d ss R ^ r 4f 
" q e Jr 

dip 7 



(B.H) 
(B.12) 



dv 2 <91ogA 

2ilo— + c — 

or or 



'dlogM 2 



dr 



v dr 
- (1 + M 



+ 

v 



The radial part of the Euler equation, together with the con- 
tinuity equation, lead to the following differential system: 



df iLoM 2 f 

v— H — 

dr 1 — M 2 

dg iuoM 2 g 

v— H — 

dr 1 - M 2 



icov 2 g . 2 SSr iu 
+ iloc e 



M 2 
iuif 



Jr v 



bo 
bi 



-iA 



c 2 (l-M 2 ) UJ 



i . , . i8K R u 
+ -Ag :V f H — e 



(B.13) 

J R v 

1 

(B.14) 



uc 2 (l - M 2 ) 

i'yA 
~«c 2 (l -M 2 ) ( 



LU + llUJ— — h L —. 

dr r 2 



JR » 

dr 

to f r dc 8 
i Jr " — 



1 



c ( ■ 0lo L 



dr cur 2 A ' 



(B.21) 

(B.22) 
(B.23) 
(B.24) 
(B.25) 



The function A is introduced by our choice of writing the dif- 
ferential equation satisfied by the pressure perturbation, which 
is justified from the physical point of view. Note that A never 
vanishes as long as the real part of the frequency is different 
from zero. 



where Ag lV is the non radial part of the Laplacian in spherical 
coordinates. The constant SKr is defined as follows: 



5K ■ 



v ■ (V X w) - c 2 A e 



SS 



We deduce from Eqs. ( B.fi| -B,7) that 5 K is conserved when J 
advected: 



( d d 1 
\di +V d-r\ 



8K = 0. 



B.2.2. A more compact mathematical formulation 
(B.15) It is convenient to define (f,g) as: 

Jr T~mI * f ; (b.26) 

f r M 2 dr 

Jr T^mZ v . (B.27) 



(B.16) 



5K can be written independently of the system of coordinates, 
as the first order perturbation of the quantity K defined by: 



K 



M v 

4-7T pV 



V x w — c V x 



v „S 
- x V- 

v 2, 7 



(B.17) 



g = e •> r » g. 

The differential system Eqs. ( |B.18| - [B.19| ) is then simpler: 

df iuvg 

~dr = 1-M 2 

dg _ if_ [ u 2 _ I?_ 

dr luv c 2 (l — M 2 ) r 2 



. c 2 SSr f r d '* ., 

■ iu ^e JR Hi-m 2 ) . 

v 7 



i8Kr ii 

+ -s — e 



, r_dr 

Jr 1.(1-/ 



(B.28) 



.(B.29) 



but the conservation of K is established only in the linear ap- 

proximation. By projecting /, g, SSr and the radial component The homogeneous differential equation (g) can be transformed 



of the curl of the vorticity (V x WR) r onto the spherical har- 
monics Yi m (8, ip), which are the eigenvectors of the Laplacian, 
we obtain: 



into a more compact form using the new variable X: 

AX 
dr 



df iuoM 2 f 

v— -\ — 

dr 1-M 2 

dg iuM 2 g 
dr 1-M 2 



iu>v 2 g . 2 SS R t u 
+ iujc e 



M 2 
iujf 



Jr " 



7 
iL 2 

c 2 (l - M 2 ) ior 2 



— 7i1 + —o— e 



(B.18) 



W 



1-M 2 ' 

1 i 2 2 \ 
-5-0 (w -ui), 



(B.30) 
(B.31) 



where uji is defined by Eq. 



where we have defined L 2 = 1(1+1) for the sake of the simplicity 
of the equations. 



M 2 



•1% 



u 8Sr I uj 



+ i 



d 2 f 
dX 2 
d\o%M 2 
dr 



+ Wf = 

SKr 



+ 



(B.32) 

With this new variable, the sonic point r s corresponds to X — > 



B.2. Second order differential equation 
B.2.1. Pressure perturbation 

The most natural function to describe the behaviour of acoustic +°°> while the spatial infinity corresponds to X = 
waves is the pressure perturbation Sp/p, which is related to the 
sound speed perturbation through: 



Sp 

IP 7 - 



2 8c SSr iw f r 
e J r 11 . 



1 c 



(B.20) 



Appendix C: Approximations of the 
homogeneous solution 
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C.l. WKB approximation jar from the accretor 

The general solut ion / of the homogeneous equation associ- 
ated to Eq. (B.32) is a linear combination of ingoing (/+) and 
outgoing (/_ ) waves, which can be approximated by the WKB 
method: 



f ~Wi exp { ±l J T^f 2 



dr 



This approximation is asymptotically valid for 

<91ogW i 

< W 2 , 



dX 



(C.l) 



(C.2) 



which is satisfied at high frequency or far from the accretor. 
The Wronskien W associated to the couple (f + ,f~) of solu- 
tions is: 



»-r'-£-r'-£ 

2iuv 



M 2 



(C.3) 
(C.4) 



The pressure perturb atio n Ap associated to a perturbation / 
is deduced from Eq. (B.2): 



7P c z 



.vi- 



/ iv df 
1-M 2 + ~udr 



(C.5) 



The normalization factor a;' in Eq. (C.l) is such that the 



press ure perturbations Ap± associated to these solutions by 
Eq. (B.2) have the following asymptotic behaviour: 



Ap ± 

7P 



Mi 



c 1±M 



(C.6) 



Both carry the same acoustic flux F — Mocto according to 
Eq. ©. ' 



C.2. Definition of the refraction coefficient 1Z S and 
estimate at high frequency 

C.2.1. 7 < 5/3 

A Frobenius expansion in the vicinity of the sonic point (M = 
1) shows that any solution of the homogeneous equation can 
be projected onto the couple ft, fa of solutions, where fi(r) is 
regular and /2(f) is singular: 



/i(r) ~ l + 0(r), 
f 2 (r) ~ [1 + C(r)]exp 



■ log(r - r B ) 



(C.7) 
(C.8) 



. 2c B M B 

This defines a unique coefficient 1Z S such that the solution 

/„ = /++£,/- (C.9) 

is regular at the sonic point. The coefficient \1Zs\ 2 is the refrac- 
tion coefficient of sound waves. The linearity of the equations 
guarantees of course that the solution Apo = Ap + + TZ s Ap~ 



is also regular. 

An expansion of W in the vicinity of the sonic point leads to 

go +qi(r - r s ) + 0(r - r s ) 2 , 
, .2 



W 

qo 



(CIO) 
(C.ll) 



'11 



Cs dr 
'2M 



+ 



2L 2 M B 



exp 



X 



(C.12) 
(C.13) 



The WKB approximation of / is valid at high frequency u 2> 



wf ut , for r - r s < 1: 



/ ~ c s exp ±i 



2c s A^f 



log(r-r s )+C(l) 



(C.14) 



The asymptotic solution to the differential equation (|9() can 
also be approximated by a Bessel function in the vicinity of 
the sonic point, for r — r s <^C 1. The regularity of /o at the 
sonic radius implies: 



/o ~ XJ V 



-c s q{ 



(r- 



(C.15) 



(C.16) 



M s CsM s 

The Bessel and WKB approximations are both valid in the re- 
gion where lo^ 1 r — r B <C 1. In this region, the matching 
of the two approximations gives the normalization constant 
A and the refraction coefficient TZ S . The asymptotic behavior 
of the Bessel function, when both the imaginary order v and 
the argument tend to infinity, is obtained from Watson (1952, 
Chap. 8.6, p. 262). The Bessel function J v (z) is written in terms 
of the Hankel functions H$p(z) and Hi> (z). Let x,z be real 
numbers, with \x\, z 3> 1, such that the order v is purely imag- 
inary v = izsinha;: 



Mz) 
H?\z) 



(C.17) 



H, 



(2), 



2 exp(7rz sinh x) 
nz cosh x 

2 exp( — 7T2 sinh x 
nz cosh x 



3 +i [z (cosh x — x sinh x) — ^-J 



(C.18) 



-i \z (cosh x — x sinh x) — ^ j 



C.19) 



Applying these asymptotic expansions to Eq. (C.15), with 



M s 



sinh x 



(r-r s ) s > 0, 



(r - r s )- 



(C.20) 



dr 



<0, 



(C.21) 



we obtain A and l/o^s)! to first order at high frequency 
27TCsO; 



Ms 



(C.22) 



|/o(r.)| = cs- (C23) 
The expon ential decrea se of \7t s \ at high frequency is deduced 
from Eqs. (C.18) and (C.l£), by identifying the Hankel func- 
tions with the ingoing and outgoing acoustic waves: 



\TZ S \ <x e c 



(C.24) 



12 



T. Foglizzo: What does perturbed Bondi accretion sound like ? 



C.2.2. 7 = 5/3, I = 

The WKB approximation is valid in the range r 3> 1/uj. 



f ~ cM 2 exp ±itu 



M 



o v(l-M 2 ) 



dr, 



(C.25) 



where the normalization is chosen in agreement with Eq. ( |C.l| ). 
In the region r < 1/lo where the WKB approximation is not 
valid, we approximate the solution by the Bessel function Jo, 
which is regular at the origin: 



&f± + - 
dr 2 r 



1 + 0(7-3) 



dfo , cu 2 ~ 



l + 0(r 2 ' 



fo ~ /o(0)J 



2- 



0, (C.26) 
(C.27) 



The other Bessel solution, Yq, is singular with a logarithmic 
divergence, as in the case of 7 < 5/3. The Bessel approximation 
for r 3> 1/u) matches the WKB approximation: 



/o(0)Jo ( ^ 

\22 



/o(0) ( — 

1 1YLUT 



2" 



2 2 2 

From this we deduce that: 



e 22 +TZ s e 22 



/„(0) = ( ' 
= 1 



(C.28) 
(C.29) 

(C.30) 
(C.31) 



C.2.3. 7 = 5/3, I > 1 

The case of non radial perturbations needs to be treated sepa- 
rately, because the behaviour of the solution is quite different 
near the sonic point. If I > 1, there is always a turning point 
beyond which the solution is evanescent, such that it decays ex- 
ponentially on the sonic point. The leading behaviour of fo, go 
on this essential singularity is the following, within a multi- 
plicative constant: 



fo 
9o 



I 2L 
r s exp r , 

T 4 

3 iL 2L 
22 — r exp -. 

lots r& 



(C.32) 
(C.33) 



The solutions of the other branch diverge exponentially on the 
sonic point: as in the case of a regular sonic point, there is a 
unique combination of ingoing and outgoing waves such that 
the solution does not diverge at the sonic point. The t urnin g 
point r t corresponds to the radius where W = in Eq. ( B.32j ) : 



,.,2 



(C.34) 



The WKB approximation is valid for r t « r « 1: 



f 



[1-££(1-A*»)]' 
M 



exp ±iu 



v(l-M 2 ) 



1-^(1-X 2 ) 



d<C35) 



Appendix D: Acoustic efficiencies Qs, Qk of 

entropy and vorticity perturbations 

D.l. Acoustic efficiency Qk 

The outgoing pressure perturbation produced by the advec - 
tion of a perturbation 8Kr is computed by solving Eq. ( B.32 ), 
imposing the regularity at the sonic point and the absence of 
an i ncomi ng sound wave from infinity. The general solution of 
Eq. (B.32) with 5Sr = can be written as: 



/ 



5p 

IP 



iSKh 
2lo 



f + a + 



r <p+ 



}+ to r _^ 



r j- iw r *r 



iSK E 

277 

IP 



r 

Ap~ 

IP 



a + 



f iu r — : 



M 2 ) 



(D.l) 



(D.2) 



Accordi ng t o Eqs. (C.C) and (C.7-C.8), the integrals involved 
in Eq. (D.2) are well defined both at infinity and at the sonic 



radius. The condition that no sound wave comes from infinity 
leads to: 



(1-A1 2 ) _ 

The regularity at the sonic radius can be written: 
(3 = -n s a- / 'A^r^I^T 



(D.3) 



(D.4) 



Thus the efficiency of sound wave emission by the advection of 
vorticity perturbations is: 



Qk 



1 

2^ 



JR 1-M 2 v 



dr. 



(D.5) 



This integral converges as 1/r 3 at infinity. It can be trans- 
formed using integrations by parts together with the homoge- 
neous differential system satisfied by fo,go, so that it converges 
rapidly at infinity, thus allowing an easier numerical computa- 
tion. The following expression converges as 1/r 5 : 



1 r tofi 



+ M- 1 dr 
1 ~V~ 



d ( 1 - M 2 \ . iuv . , 
^1 J /o + — go * dr - (D ' 6) 



D.2. Acoustic efficiency Qs 

The outgoing pressure perturbation produced by an incom- 
ing entrop y perturbation SSr is also computed by solving 
Eq. (B.32), imposing the regularity at the sonic point and the 
absence of an in comin g sound wave from infinity. The general 
solution of Eq. (B.32) with SKr = can be written as: 



- _ SSr 
2 7 
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fl-M 2 tof-jJt 



dr V A4 2 



(D.7) 



Sp _ SS R 
7P 27 



Ap~ 
7P 
Ap+ 
7P 



P+ f 



9 ( i-M 2 r r^, 



dr \ M 1 



a+ f 



Or \ M 2 



{l-M 2 ) 



(D.. 



Since the integrals of Eq. (D.J) converge at the sonic point, the 
regularity of the solution at the sonic point can be written: 



Two successive integrations by parts make the convergence of 
the integral much faster, thus allowing an easier numerical cal- 
culation: 



-4. 



o_ 

dr 



+ 



c£\ d 



W / or 



(c -V ) 



35 = Tr< {1 ~ M) 



v d ^ 
icu dr 



V* 1- 



, (D.18) 

D.19) 



(3 = -alls 



fa- 



d ( 1-M\ i»r^ 



dr \ M 2 



(D.9) 



Appendix E: Asymptotic estimate of |Qs|, \ Qk\ 
at high frequency for 7 = 5/3 



In Eq. ( D.8 ) the integrals diverge at infinity. Two integrations 
by parts leads to define the two functions F±,H± as follows: 



Eq.(D.14) can be estimated in the limit ui 3> c^/GM, when 
7 = 5/3. We treat the case of radial (I = 0) and non radial 
(I > 1) perturbations separately. 



M 2 



IP 



(D.10) 



dr 



v(l-M 2 ) 2 dp 



iioM 2 Or 



E.l. Radial perturbations I = 
(D.ll) w e use a local expansion of M where r«l: 



The pressure perturbation associated to the entropy perturba- M — 1 — 4r 2 + 0(r). 
tion SSr is: 



(E.l) 



5p 

IP 



+ 



SS R 

27 

Ap+ 

7P 



7P 



F + (r Q )+ / H + dr 



Writing /o = /- + lZ 3 f+, we obtain for ri > and r2 < 1: 
rr2 J«f:^^ d (v{l-M 2 ) 2 0~f Q 



F~(r )+ / H'dr 



(D.12) 



Or 



M 2 



dr 



6* 



Far from the accretor, the pressure perturbation can be pro- 
jected onto the ingoing and outgoing sound waves, a is chosen 
such that there is no sound wave coming from infinity: 



dr 



— du, 



I 3 
2^r2u 



a = F (r ) + / H'dr. 

'To 



(E.2) 



,1/2 



(D.13) 



The range of validity 1/uj <C r <C 1 implies that u\ S> l/i 
and ti2 < The contribution of the region r^ < r < 00 is 
negligible. 



Using Eqs. (D.9) and (D.13), and the fact t hat F (r 3 ) = 0, we 
take the limit r — » 00, ro — ► r s in Eq. (D.12) and obtain: 

Qs 



= / e J i,(i-.m 2 ) — 



9r 



.M 2 <9r 



dr<aj3. (E.3) 



27LJ 



9r 



w(l - M 2 ) 2 dfo 



M 2 



Or 



dr.(D.14) contribution of the region r < n is estimated using the 

asymptotic behaviour of the Bessel function Ji : 



This integral can be written with the foll owing form at, using 
the homogeneous differential system (Eqs. B.28 , B.29 ) satisfied 
by fo,9o- 



d_ { l-M 2 df \ 2 ii /o(0) 

iA dz 



Qs = TT , 

2 7 J r , 



, f r l + M- 1 d 

'J T~f^~ 



» (Afc/o + B feS ro)dr, 



(D.15) 



9z V -M 2 dz 

d ( l-M 2 Of 



ziJi(z) 



Oz 



M 2 Oz 



. /o(0)24 _.- 
dz ~ - — rp- — r e 4 Zi 



where the functions ^4fc(r), Bk(r) depend only on the unper- 
turbed flow, and the integral converges as l/r fe far from the 
accretor: 



Thus we deduce the acoustic efficiency at high frequency: 



Ax 
Si 



ioj 

V 

0_ 

Or 



u 2 

I 2 2x 
(C -V ). 



-du. 



(E.4) 
(E.5) 

(E.6) 



« 3 



(D.16) 



The main contribution to the integral co mes from the region 
where u ~ 1, i.e. r e fi ~ (6/lo) 2 ^ 3 /2. Eq. (E.6) is rewritten in 
(D.17) Eq. (Eq) using the Gamma function. 
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E.2. Non radial perturbations I > 1 

As in the case of radial perturbations, the main contribution 
to the integral comes from the region accessible to the WKB 
approximatio n. Th e outgoing wave is the main contributor to 
the integral (D.14). Usi ng t he asymptotic development of v,c 
for r -C 1 in Appendix. A. 2, 



[\ 


1-M 


'fa 





L 2 c 2 



(l-M 2 ) 



dr 



_22 r 2 w _|_ 



v(l -M 2 ) 
L 2 



22 7-2 ^ 



(E.7) 
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With u = 2 2 r2oj/3, the acou sti c eff ici encies at hi gh fr equency 
are deduced from Eqs. flC.35j), (Euj]), flD.14p and (JeTtI): 



2s 
Qk 



7 U/ 



7 

.2*H 



exp - 



9u 



du 

U3 



exp - 



Using the new variable 



9w 



du 

5" * 

a 3 



t = 



2L 



9u 



(E.8) 
(E.9) 

(E.10) 



the integrals in Eqs. (E.8-E.E) can be written using K2, the 



modified Bessel function of order 2/3 (Gradshteyn & Ryzhik 
1980, p. 430): 



exp - 



Thus we obtain: 



L 2 \ 


du 


9u j 








9u j 


M3 



,(E.ll) 



37 



loL 2 



-K2 



(¥)■ 

^~!*.(i,) 4 ^(f 



(E.13) 
(E.14) 



Here again, the main contribution to the integrals comes from 
the region r e g ~ (6/lo) 2 ^ 3 /2. By contrast with \Qs\, \Qk\ de- 
creases at high frequency. 
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